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We examine electrolyte systems confined between two parallel, grounded metal plates using a 
field-theoretic approach truncated at one-loop order. For symmetric electrolytes, the system is 
locally neutral, and the electric potential between the plates is everywhere zero, in agreement with 
the predictions of the Poisson-Boltzmann equation. The density distribution of the ions, however, 
is non-uniform, with a maximum value at the center of the system and that lies below the bulk 
electrolyte density. This differs from the Poisson-Boltzmann result, which is uniform and equal to 
the bulk electrolyte density. For asymmetric electrolytes, the system is no longer locally neutral, 
and there is a non-zero static electric potential. In addition, the surface of the metal plates becomes 
charged. When the plate separation is smaller than the screening length, the charge distribution 
has a single peak at the center of the system. As the separation becomes greater than the screening 
length, the charge distribution possesses two peaks with a local minimum at the center of the system. 
These peaks remain localized near the surface of the plates, while the bulk of the system has nearly 
zero charge. When the plates are immersed in a bulk electrolyte solution, an attractive fluctation 
induced force develops between them that is analogous to the Casimir force. Unlike the Casimir 
force, however, this attractive interaction has a finite limit in the zero separation limit and decays 
inversely with the separation at large separations. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The electrical double layer near a solid/electrolyte in- 
terface occurs in a wide variety of contexts and plays a 
crucial role in many colloidal, electrochemical, and bi- 
ological systems. Much of our understanding of these 
systems are based on the Poisson-Boltzmann equation, 
which has been quite successful 0, 0, E|- This equa- 
tion is based on a mean-field type approximation, and, 
consequently, neglects fluctuation effects which result in 
ion- ion correlations. In many cases, these fluctuations 
are not crucial, and only add a qualitative change in the 
predictions. 

In certain situations, however, the contribution of fluc- 
tuations are important and lead to qualitatively different 
behavior than that predicted by the Poisson-Boltzmann 
equation PI . Examples include the presence of an attrac- 
tive interaction between two like-charged colloids |^ 0| , 
the collapse of polyelectrolytes 0] , and charge segrega- 
tion in confined electrolytes 

Due to the importance of fluctuations, many the- 
oretical methods have been developed to incorporate 
fluctuation corrections in electrolyte systems. Some 
of these include liquid-state theory approaches based 
on the Ornstcin-Zernike equation coupled with an ap- 
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proximation closure, such as the hypernetted-chain or 
mean-spherical approximation 0, ^i Qi density 
functional theories [ isl ll^ . the modified Poisson- 
Boltzmann equation |l7| . computer simulation meth- 
ods llEl llSL IT^ , and field theoretic approaches 

We examine a pair of grounded, metal plates which are 
separated by a distance L and immersed in an electrolyte 
solution which is dissolved in a continuum with uniform 
dielectric constant e. In this situation, the Poisson- 
Boltzmann equation predicts that the system behaves 
trivially: the ion concentration between the plates is uni- 
form, and the electric potential is everywhere zero. For 
this system, the fluctuation corrections are non-trivial 
and lead to some interesting phenomena. In this case, 
the ion concentration is non-uniform, and in some sit- 
uations can lead to the presence of a non-zero electric 
potential and a surface charge on the metal plates. In 
addition, we flnd a fluctuation induced attractive force, 
analogous to the Casimir force, between metal plates. 

The remainder of the paper is organized as follows. In 
Section m we quickly review the functional integral for- 
mulation of electrolyte systems. Then in Section ITTll we 
quickly summarize the main results in mean-fleld theory 
(i.e., the Poisson-Boltzmann equation) for electrolytes 
confined between two metal plates and then present field 
theory calculations to order one-loop. Finally, in Section 
II VI we summarize the major finding in this work. 
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II. FUNCTIONAL INTEGRAL FORMULATION 

Consider an electrolyte system that is immersed in a 
continuous medium with a uniform dielectric constant e 
and is composed of point charges. The charge density 
Q(r) in the system is given by 



g(r) - ^ <z„<5(r - r„,fc) + S(r) 



(1) 



where r^^k is the position of the fcth particle of type a, 
and S(r) is an applied charge density. 

The total electrostatic energy E of the system is given 
by 11 



E ^ — I drV0(r) • W<l>{r) 



(2) 



where (j) is the electrostatic potential, which can be ob- 
tained by solving Poisson's equation 



4^Q(r) 



(3) 



where (f> is the electrostatic potential, and Q{r) is the 
charge density. 

The associated Green's function to this problem is: 



47r 



V2Go(r,r') 5(r-r') 



(4) 



where G'o(r,r') is the Green's function. Given the 
Green's function, the electrostatic potential can be deter- 
mine directly from the charge density distribution Q(r) 
in the system j:24j : 



dr'Go(r',r)Q(r') 

+ ^ ^dr'[Go(r',r)V>(r') 
-<^(r')V'Go(r',r)] (5) 



where the surface integral is taken over the boundary of 
the system, with the vector dr' pointing outwards from 
the system, normal to its surface. The second term de- 
pends only on the boundary conditions of the system and 
not on the charge distribution within the system. 

For the problem of an electrolyte confined between two 
metal plates, we have Dirichlet boundary conditions and 
the potential 4> is zero on the metal surface. In this situ- 
ation, we choose Go(r',r) = for r' on the surface, and 
so the surface terms in Eq. vanish to yield: 



dr'Go(r,r')g(r') 



(6) 



The energy of the electrostatic field can then be written 
entirely in terms of the charge distribution 



E dvdv'Q{T)Go{vy)Q{T') 



(7) 



The expression for the energy given above contains con- 
tribution from the interaction of each ion with itself — 
the self-energy. Subtracting this term from the energy, 
we write 

^ ^drd/Q(r)Go(r,r')0(r') 



E = 



(8) 



a.k 



The grand partition function Zq for a system of point 
charges at fixed chemical potential //q, is given by 



oo 



Zg[i] 



E E n 



1 



Ni=Q Nm=0 V 

X / I I avrte 



3N^ 



Tt 



E E n 



1 



J<irdr'Q(r)Go(r,r')Q(r') 

vpEc.fc[7c(r„,fc) + |9^Go(r„.,,r„,fc)] 



(9) 



where (3 — l/{kBT), ks is the Boltzmann constant, 
is the number of ions of type v, Vaj is the position of the 
jth ion of type a, Aj/ is the thermal wavelength of an ion 
of type ly, 7a (r) = /3[/ia + Ua(r)], and Ua is an external 
potential acting on ions of type a. 

Introducing the Hubbard-Stratonovich transforma- 
tion, the grand partition function can be expressed in 

terms of a functional integral USUI mm 



Mo 



X exp 



xe2 



^ j drdr'V(r)Go-i(r,r')^(r') 

(10) 



where 

AAo = 



'Dip{-) exp 



j drdrV(r)Go-i(r,r')V(r(:iL) 



The function itp(r) can be interpreted as being equal to 
an instantaneous value of f34>{r), and the functional in- 
tegral can be thought of as an integral over all possible 
"shapes" of the electric field due to the thermal motion 
of the electrolytes. 



A. Mean-field approximation 

In the mean-field approximation, the value of the func- 
tional integral is replaced by the value of the integrand 
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at its saddle point. The value of the field at the saddle 
point, denoted by "0, is given implicitly by 



5 In Zg [iip] 



0, 



which leads to the relation 

a 

+E(r). (12) 



This is precisely the Poisson-Boltzmann equation. The 
electric potential is given by (j){r) = kBTiip{r). 

Within the mean-field approximation, the grand po- 
tential is given by 

InZcb] = ^/ drdr'z^(r)Goi(r,r')*^(r') 

+ ^K'^ f dre''°('')-'°*'^('") (13) 

The density distribution of ion of species a can be deter- 
mined by taking the functional derivative of the grand 
partition function with respect to 7^ j^], which yields 



S In Zg [7] 
'57a (r) 

^-<ig7c«(r)-gQii/'(r) 



(14) 



where the n-th moments of the mean-field charge densi- 
ties have been defined as 



Kn{r) 



(16) 



The second moment if 2(1") corresponds to a (in general 
r-dependent) inverse screening length. The third mo- 
ment Ksij) will turn out to be crucial for the physics of 
asymmetric electrolytes. Expanding to order one-loop, 



In [7] = In^cM 
1, 



(17) 



B. One-loop expansion 



- -Tr[ln(l + if2Go)-Mo]. 

Here, the last term iTr[^2Go] corresponds to the self- 
energy of the charges which has to be subtracted from 
the total electrostatic energy. Equation p7|l . which is 
still completely general, forms the starting point for the 
discussion of fluctuation- induced effects in electrolytes. 
Since in the following, all the calculations are to one-loop 
order, the index (1 ) in Z^^^ is omitted. 



III. CONDUCTING PARALLEL PLATES 

Consider a medium of dielectric constant e which is 
confined between two conducting metal plates separated 
by a distance L. In this geometry, the mean-field solution 
for the field ■0(r) and the densities Pa(r) is trivial. 



The functional integral for the grand partition func- 
tion can be cast into a form that allows one to perform a 
loop- wise expansion around the mean-field solution ipi^). 
Defining 5ip = iJj — ip a.s the deviation from the mean-field 
solution, Eq. is re-written as a functional integral 
over dip(r). The quadratic term in the action is now de- 
fined by a screened propagator (Green' s function) defined 
by a Dyson equation 

G-i(r,r') - G^\r,r')+5{v~r')K2ir'), (15) 



i^ir) = 0, p,(r) =p„ = A-'^eT" 



(18) 



That is, the electric potential between the plates is zero, 
and the electrolyte density is uniform and solely deter- 
mined by the chemical potentials /Zq, and the temperature 

ksT = 1//3 (7q = /3fia). 

The (screened) Green's function G(r, r') for the 
Poisson-Boltzmann equation is obtained from a Fourier 
series for the screened Coulomb interaction and given by 



G(r,r') = ^ly 

n—l 



1l 



fn'KZ\ 
sm(— jsm 



Jq-L(x-x')+iq2{y~y') 



(19) 



where n is the inverse Debye screening length in the 
mean-field approximation, which is defined as 



47r 



47r/3 



YllPa 



(20) 



The (unscreened) Green's function Go(r, r') of the Pois- 
son equation (see Eqs. Q and Q), is obtained from 
Eq. I^ni) by setting k = 0. 
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A. Density profile 

All the thermodynamic properties of the system can be 
derived from the grand potential. The density profile to 
one loop order is given by the functional derivative 



Pair) = 



d In Zq 
<57a(r) 



Pa 



1 



AG(r,r) 



(21) 



2 



dr'AG(r',r')i^3(r')G(r',r) 

^ J 

where 

AG(r,r) = G(r, r) - Go(r, r 
with Ii defined as 



—h(KL,z) (22) 



Ii{kL,z) 



^sin^ 



mrz 

IT 



In 1 




Note that in the bulk limit where L oo, G(r,r') = 
g-«|r-r I ^ _ g^j^j AG(r, r) = —n/e. In this case, 



the density profile becomes uniform and is given by 



2e 



2eK 



(24) 



For symmetric electrolytes , where the magnitude of all 
the charges of the particles is the same (i.e., \qa\ = q for 
all a, and therefore = 0), the last term in Eq. H22fl 
vanishes, and the expression for the density profile re- 
duces to 



1 



eL 



Ii{kL, z) 



(25) 



The second term is the only correction term for a sym- 
metric electrolyte, and its sign is independent of the sign 
of the charge. As a result, the fluctuation corrections 
increase the electrolyte concentration of the system, as 
compared to the mean-field approximation. However, 
the ion concentration at the metal surface remains at its 
mean-field value. In addition, the system is everywhere 
locally neutral. 

The difference between the density between the plates 
and that of a bulk system at the same chemical potential 
is given by 



Pa(r) -p, 



bulk 



eL 



Ii[kL,z) - — 



(26) 



Plots of the density profile difference between the con- 
fined system and the bulk system at different values of 
kL are given in Fig.Q] In Fig.^, the plate separation L 
is varied while holding the value of k (i.e., the ion chem- 
ical potentials) fixed; in Fig. the value of k is varied 
while the plate separation is held fixed. 

The ions are more concentrated in the center of the 
system than near the plates. At the center of the plates, 
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FIG. 1: Density profile for a symmetric electrolyte solu- 
tion confined between two conducting plates with (a) k held 
constant {pa ~ fipaq^i-^/ €) and with (b) L held constant 
(pb = Ppaqa/i^L)): (i) kL = 0.1 (solid line), (ii) kL — 1 
(dashed line), (iii) kL = 10 (dotted line), and (iv) nL — 100 
(dashed-dotted fine). 



where the ion concentration is at a maximum, the density 
can be computed in closed form, to yield: 



Pa 



1 



In cosh — 
eL 2 



(27) 



In addition, for a symmetric electrolyte, the concentra- 
tion of the ions in the confined system is lower than that 
of the bulk system at the same chemical potentials. 

The local charge distribution Q(r) between the plates 
can be computed from the density distribution of the 
ions: 



Qir) 



J2qaPa{r) (28) 

a 

I AG(r, r)-K2 j dr'AG(r', r')G(r', r) 

^i/l(Ki,z)-(Ki)2/2(Ki,z) 



sinh^sinhij^(z/L - 1) 
cosh^^ 



where we have used the fact that J2a la Pa ~ 0, de- 
fined cr* = K^/{2e), and introduced the fimctions (cf. 
Appendix El 



I2{kL,z) = ^ 



In 



(i + (^y 

(2m7r)2 + {nLf 



(29) 



{nLy 



E 



(2m.)Mn(l+(^)^ 
(2m7r)2 -t- (kL)2 



In the case where = (e.g., a symmetric electrolyte 
with \qa\ = q for all a), Q(r) is trivially zero, and the 
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FIG. 2: Variation of the charge density for an electrolyte 
solution confined between two conducting plates at constant 
k: (i) nL = 0.1 (solid line), (ii) kL — 1 (dashed line), (iii) 
kL — 10 (dotted line), and (iv) kL = 100 (dashed-dotted 
line) . 



system (apart from being locally neutral everywhere) is 
globally neutral. 

A non-trivial situation arises for asymmetric elec- 
trolytes with K3 0. As shown below, the fluctuations 
around the mean-field solution lead to a non-zero charge 
distribution between the metal plates with a density pro- 
file that strongly depends on the dimensionless parameter 
kL. 

In Fig.|51 the variation of the local charge distribution 
for different values of kL is plotted at constant k. For 
low values of kL, the local charge distribution has a single 
peak located in the center of the system with the same 
sign as K3. For a binary asymmetric electrolyte, this 
implies that the sign of the net local charge is the same 
as that of the ion with the larger charge magnitude. 

As kL further increases, the charge distribution ini- 
tially uniformly increases; however, at a critical value 
of kL, the single peak becomes a local minimum, and 
two peaks emerge on either side of this minimum. These 
peaks move towards the boundaries with increasing kL, 
while the value of the charge density at the center of 
the system decreases, eventually changing sign. This is 
similar to the change in the monomer distribution for a 
polyelectrolyte confined between two oppositely charged 
plates J7]. 

As the plates become further and further separated, 
the charge density near the plates approaches a limit- 
ing form Qoo(z) = lim.KL^oo Q{z)/(2Ka*), which can be 
written analytically as: 



Qoo{z) 



1 



-2kz 



1 logs 



[Ei(-3Kz) - Ei(-Kz)] 

-.z 

- [Ei(-Kz) - Ei(Kz)] 



(30) 




FIG. 3: The local charge distribution near a metal plate in 
the limit of infinite plate separation. 



where Ei denotes the exponential integral. This function 
is plotted in Fig. O 

Integrating the distribution of charge contained be- 
tween the plates, the total net charge per unit area Q 
is determined to be 



Q = 



/ dzg(r) = 2cr*(Ki)tanh — h{nV). (31) 
Jo 2 



Therefore, an asymmetric electrolyte system has a non- 
vanishing net charge. As shown in the next section, 
this net charge is counter-balanced by an induced sur- 
face charge on the metal plates. 



B. Electric potential 

The electric potential (/)(r) between the plates can be 
obtained from the grand partition function [2^ 



2 

87rcr*L 



dr' l\G{r' ,r')G{r' ,r) 



sinh ^ sinh^(z/i - 1) 
cosh^ 



-l2{nL,z) 



(32) 



In the case of a symmetric electrolyte, the correction term 
vanishes, and the electric field is identically zero between 
the plates, similar to the mean-field approximation. In 
the case of an asymmetric electrolyte, however, the fluc- 
tuation correction is non-zero, and the electric potential 
varies within the cavity. 

In Fig. ^ the electric potential 4>{z) in the system is 
plotted. If the plate distance is kept fixed, then the elec- 
tric potential between the plates uniformly increases with 
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FIG. 4: Variation of the electric potential between the metal 
plates at (a) fixed k (0a = —8TTa*/{eit)) and (b) fixed plate 
separation L {(j>h — —8na*L/e): (i) kL = O.l (solid line), (ii) 
nL — 1 (dashed line), (ill) kL — 10 (dotted line), and (iv) 
nL = lOO (dashed-dotted line). 



K to a constant profile (see Fig. 2Ji). If k is field fixed 
while the distance between the plates L is increased, then 
the potential increases from zero to a maximum value at 
a critical value of L, and subsequently decrease to zero 
(see Fig.UlD). 

As a result of the different density distribution between 
ion species, local charge densities develop. The elec- 
trolyte system is no longer necessarily electrically neu- 
tral overall. The metal plates, however, develop a surface 
charge, which keeps the overall system globally neutral. 
The magnitude of this surface charge is given by 

e 90(r) 



z 



z=0 

a*{nL) tanh ——I^{kL) 



(33) 



Note that the surface charge density on each of the plates 
is precisely half of the net charge of the electrolytes be- 
tween the plates (i.e., a = — Q/2), so the system as a 
whole is electrically neutral. 

In Fig. |S1 we plot the variation of the surface charge 
on the metal plates as a function of kL. As kL increases, 
the surface charge increases, reaching a constant non- 
zero value a/a* = I — ln3/2 sa 0.45i as kL approaches 
infinity. This implies that in the limit where the plates 
arc infinitely far apart, the fluctuations in an asymmetric 
electrolyte still induce a surface charge. 



C. Casimir Force 

We now consider the situation where two, grounded 
metal plates are immersed in a bulk electrolyte solution. 
The chemical potentials of the electrolytes in the bulk 
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FIG. 5: Variation of the surface charge density with kL. Note 
a* = -K3/[2e). 



solution are fixed at 7^ = /J/Iq, which are identical to 
the chemical potentials of the electrolytes between the 
plates. We assume that the thickness of the metal plates 
is much larger than the size of the electrolytes and much 
larger than the separation between the plates. 

In order to obtain the force between the two plates, we 
calculate the pressure which is given by 



P 



dlnZa kBTdlnZa 



dV 



A dL 



(34) 



where A is the surface area of the plates, and V = AL 
is the total volume of the system. We therefore need the 
full expression for the grand potential Zq, which to order 
one-loop is explicitely obtained from Eq. (|i7|l and can be 
written as (see Appendix IB|) 



InZcM = AlY^Pa 

a 

2 



"8^ 



■I{kL) 



I{kL) = - 



( sinh; 



dxx In . 
V ^ 



(35) 



For large values of L, the the grand potential behaves as 



AL—-A\ ln2- - 
I27r V 2 



2 

8^ 



^2 , A C(3) , , 

-A-ln(.L)--iU + ... (36) 

where the higher order terms vanish more rapidly than 
with increasing L. 
The first two terms are proportional to the volume of 
the system. This is essentially the Debye-Hiickel approx- 
imation (order one-loop) for the fluctuation correction 
for a bulk electrolyte. The third term is proportional to 
the surface area of the metal plates and is independent 
of the separation distance L. This term is esentially the 
free energy of creating the metal/electrolyte interfaces. 
The higher order terms represent the free energy cost of 
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FIG. 6: Fluctuation induced attraction between two parallel 
metal plates separated by a distance L. 



separating the plates. These terms are responsible for 
inducing a net attraction between the metal plates, as 
demonstrated next. 

Owing to Eqs. H34|l and (|35|) . the pressm'e in between 
the metal plates, due to the presence of the electrolytes, 
is given by 



^3 

Y^p^-—I'{nL 



(37) 



Using the large-i expansion Eq. H36|) . for kL ^ 1 the 
pressure behaves as 



C(3) 1 

12^ ~ 8^L ^ "8^l3 



(38) 



The sum of the first two terms on the right-hand side 
are equal to the pressure of the bulk solution (i.e., the 
solution outside the plates). This pressure acts on the 
outside surface of the plates. 

The net attractive force per unit area V acting between 
the plates is given by the difference in the pressures inside 
and outside the plates: 



ViL) = piL,j) -p(oo,7) 



Sir 



I'{kL) + - 



(39) 



In Fig. IHl the variation of this attractive force with plate 
separation is plotted, while holding the inverse screening 
length K constant. 

For small values of kL (e.g., at small plate separations), 
the force between the plates behaves as 



V{L) 



127r 



1 - 



kL {nLf 



90 



(40) 



The leading term is proportional to the square root of 
the temperature. In the limit that the plate separation 



vanishes, we see the force approaches its maximum value, 
which is finite. 

For large values of kL (e.g., for large plate separations), 
the force behaves as 



.2 



V{L) 



1 - 



C(3) 

(kL)2 



(41) 



The leading order term is independent of temperature. 
For large plate separations, the force decays inversely 
with the plate separation. 



IV. CONCLUSIONS 

We have performed order one-loop calculations for an 
electrolyte confined between a pair of grounded, metal 
plates. For this system, we find that the fluctuation cor- 
rections lead to qualitatively different physics than pre- 
dicted by the mean-field approximation. 

The density profiles of the ions are non-uniform. For 
symmetric electrolytes, the densities of the ions have 
their minimum values (give by the mean-field values) at 
the metal interface and are peaked in the center of the 
system. The magnitude of the density profiles is always 
lower than the densities of the ions in an equilibrium bulk 
solution. Despite the non-uniform density distribution, 
the system still obeys local electrical neutrality, and the 
electric potential is everywhere zero. 

For asymmetric electrolytes, local charge separation 
can occur. When the plate separation is much less than 
a screening length, the charge profile has a single peak at 
the center of the system. When the plate separation be- 
comes much greater than the screening length, the charge 
profile develops two peaks, with a minimum in the cen- 
ter of the system. This non-uniform charge distribution 
generates a non-zero electrical potential throughout the 
system. In addition, a surface charge arises on the metal 
plates. In the limit of infinite plate separation, this sur- 
face charge approaches a finite limit a/a* « 0.451. In ad- 
dition, the charge distribution near the plates approaches 
a limiting form. 

When the plates are immersed in a bath of electrolytes, 
the ion-ion correlations induce an attractive force be- 
tween them. This force is greatest when the plates are 
touching and decays as 1/L at large plate separation. 

In these calculations, the size of the ions was neglected. 
Consequently, the predictions are valid only in the limit 
where the plate separations are much greater than the 
ion sizes. Excluded volume interactions, as well as other 
types of forces, can be incorporated into the theory by 
changing the starting partition function Eq. Q [25ll2a |. 
This method can also be used to study other types of elec- 
trolyte systems, such as polyelectrolytes or even quantum 
systems. This work is currently being pursued. 
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APPENDIX A: INTEGRALS FOR DENSITY 

For the evaluation of the charge density profile in the 
asymmetric case ^ 0, Ea. H28|l . one requires the fol- 
lowing integrals: 



/ dz'AG{z',z') sin 

Jo L 



(Al) 



m — 1 



(2m7r)^(l — cosnTr) 



mTT J I n7r[(n7r)^ — (2m7r)] 



and 



J dr'AG(r',r')G(r',r) 



(A2) 



sm ■ 



&ttL 



dz'AG{z', z') sin ■ 



^ (2m7r)2 In f 



m—1 

oo 



kL 



(f — cosnTr) sin 



^ 7i7r[(n7r)^ + (kX)^] [(riTr)^ — (2rn7r)] 



To perform the summation over the index n, we first note 
the following relation, 



E 



n7r(l — cosnTr) sin 



riTTZ 

L 



(n7r)2 + a? 



cosh§(2z/i - 1 
2 cosh(a/2) 



U3) 



which can be verified by expanding the right-side of the 
expression in a Fourier sine series. In addition, we note 
that 



(A4) 



(n7r)[(n7r)2 + (kL)2] [(n7r)2 - (2m7r)2 



(2m7r)2 -f (kL)2 \ {nLf 
I 



(2TO7r)2 



TITT I 

(n7r)2 -|- (kL)2 rivr 
1 



(n7r)2 — (2m7r)2 titt 



This is used when performing the summation. 



/ 



(2m7r)2 
2 ^ (2m7r)2 + (kL)2 

m— 1 



dr'AG(r',r')G(r',r) = ^ 




cosh^^(2z/i- 1) 



_J 

cosh^ 



2m7rz 
cos 1 



e2 ULY ^ (2?7i7r 

^ ^ m—1 ^ 



sinh sinh ^(z/i - I) 
' cosh^ 



V 2m7r 



miTZ 



SvrL fsinh^f sinhi^(z/L-I) ^ (2"^'^)' In ( 1 + (^)^ 



-E 



(Ki)2 cosh ^ 
Infl-K(^)' 

\ V m-TT / 



E 



(2TO7r)2 -h (kL)2 



SttL 



' (2TO7r)2 + (Ki)2 L 

sinh ^ smh^{z/L - 1) 



cosh^ 



/3(kL) - hinL^z) 



(A5) 



where /2 and I3 are defined in Eq. H29|l . 



Finally, for the electric potential, we require the deriva- 
tive 



(A6) 



AttkL kL 
5— tanh— /3(kL) 
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which can easily be obtained by differentiating Eq. I|A5|I . 

In order to perform asyptotic analysis for the functions 
Ii, I2, and /s, it is convenient to express them in terms 
of integrals rather than an infinite series. This can be 
acheived by first noting that 



ln|l+(^ 

nir 



and then using the relation: 



2xdx 



(n7r)2 + 



xdx 



1 



sinhi^(l - z/L)sinh^^ 



sinh^^ 



sinh^^(l - z/L) sinh 



KXZ 

2 



X sinh ^ 



(A7) hi^L) = -1- 



xdx 



'2kL Jq x'^ — 1 

, kLx kL 

X x coth — coth — — 

2 2 



(AlO) 



(All) 



cos2m7ra; cosh ^(1 — 2a;)Q; 1 

4asinh(Q!/2) 2a'^ 



(2m7r)2 + 
It can then be shown that 



(A8) 



Ii{kL, z) 



{kL) /"^ sinh ^^(1 - z/L) sinh 



dx- 



sinh ^ 



fA9) 



APPENDIX B: ONE-LOOP CORRECTION 

The one-loop correction to the grand potential is given 
by the expression 



SlnZch] = --Tt[\ii{1 + K2Go)- K2G0] 
f dqi dq2 ( ( 



[n^/LY +ql + ql) [n^ / + qf + q^ 
I 



I(kL). (B1) 



We need to evaluate the following integral 

00 ^cx; 



I{kL) ^ ^ dx 



In 1 



1 



{nil / kLY + X 



(tivt/kL)^ + x 



(B2) 



n=l 



In 1 



In 



KkL) 



X 



1 



(B3) 



Now, use the relation 



da;ln(l + x/a^) = (1 + a^)ln(l + 1/a^) - 1 (B4) 



/o 
to find 



f ( 



x{KLf 



V \/xkL 



W Jo 



dxx In 



/ sinh a; 



(B5) 



Asymptotic limits can easily be obtained from this ex- 
pression. In the limit kL — > 0, we find 



I{kL) 



24 1080 



12 



540 



and 



(B6) 



(B7) 



On the other hand, in the limit nL — > cxd, we find 



I{kL) = -^(kL) + nn2- i +ln(KL) ) (B8) 



C(3) 



and 



I'{kL) 



{kLY 



2 1 C(3) 

3 kL {kLY 



(B9) 
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